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ABSTRACT 

Supernovae remnants are known to be luminous sources of infrared [Fell] emission. By 
studying how the luminosity scales with age, environment and other relevant factors, we 
can construct an [Fell] luminosity function for supernovae remnants. This will enable 
us to predict supernovae rates in starburst galaxies that are too distant for individual 
remnants to be resolved. First, however, we require accurate luminosities for a sample of 
remnants of varying ages, and in varying physical environments. As part of this project 
we have carried out an initial study of a small sample of evolved (ages greater than a 
few thousand years) remnants in M33. From these data we tentatively conclude that 
there is evidence for the peak luminosity in the [Fell] lines of these sources to arise in a 
narrow range of ages. In other respects, the M33 remnants are similar to their galactic 
and Magellanic Cloud counterparts in the observed peak luminosity. From this, and 
internal evidence as to the environment present in these regions, we conclude that the 
luminosity of evolved remnants is only marginally dependent on density and metallicity. 



INTRODUCTION 



Q | t has been realised for some time now that the intense [Fell] 
I emission observed in the near infrared in starburst galaxies 

P is likely to be due to supernovae remnants (SNR). For ex- 
sample, Moorwood & Oliva (1988) detected 1.644/tm [Fell] 

^ emission from about half their sample of galaxies classed 
as star forming. More recent work has only confirmed this: 
in their study of compact sources in M82, Greenhouse et 
al. (1991), show that the observed 1.644/tm to Br7 ratio in 
the known radio SNR, 44.0+59.5 and 41.9+58.1, is greatly 
enhanced compared to the HII regions in the same galaxy. 
Similarly, Forbes et al. (1993) find a correspondence between 
some of the non-thermal compact radio continuum sources 
in NGC253 and enhanced [Fell] emission. This outline is also 
in agreement with infrared studies of SNR in our own galaxy, 
and in the Large Magellanic Cloud (eg. Oliva, Moorwood & 
Danziger 1989). 

Some attempts have been made to use this fact to pro- 
vide constraints on the nature of the initial mass function in 
starburst galaxies, most notably the work of Colina (1993). 
This parallels similar attempts to use high resolution radio 
data for the same purpose (e.g. Colina & Perez-Olea 1992) 
under the assumption that the non-thermal emission seen 
in luminous starburst galaxies arises in SNR. However, in 
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principle, the use of [Fell] as a quantitative indicator of the 
supernova rate is superior to the radio technique, since the 
[Fell] emission from HII regions is very weak, whereas ther- 
mal radio emission from such sources can compete with non- 
thermal emission from SNR as the dominant source of radio 
luminosity in starburst systems (eg Condon et al. 1991). Ad- 
ditional uncertainty arises in the radio work since we cannot 
even be sure that the non-thermal emission seen does not 
arise in a low luminosity active nucleus (Lonsdale, Smith & 
Lonsdale 1993). However, in order to achieve our aim of us- 
ing the [Fell] lines as diagnostic probes, we require a greater 
understanding of the evolution of the [Fell] line luminosity 
as a function of age in SNR in differing environments. For 
example, the brightest compact SNR's observed in M82 have 
1.257/tm [Fell] luminosities of ~ 10 32 W, whereas a typical 
bright remnant in the LMC has a luminosity of ~ 10 30 W 
(Oliva et al. 1989). In order to derive useful supernovae rates 
from even more distant galaxies in which we cannot resolve 
the individual remnants, we must be able to predict a typical 
SNR luminosity given that we know its age and something 
about its environment. Given a supernovae rate we can then 
place constraints on the stellar initial mass function, since 
the luminous radio supernovae are typically type II (or lb), 
and hence their progenitors are massive stars with initial 
masses > 8 Mq. 

In order to calibrate this method we have therefore em- 
barked on a series of observations of individual SNR. It is 
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envisaged that the complete sample will span a wide range 
of [Fell] luminosity and environment. We require accurate 
distances in order to derive accurate luminosities. Since the 
distance to most galactic sources is poorly defined, we have 
excluded such sources from our sample. Although there are 
nearby starburst galaxies in which individual bright SNR 
have been resolved by radio observations, thus giving us 
access to the high luminosity end of our sample (eg. M82 
Huang et al. 1994, NGC253 Ulvestaad & Antonucci 1991, 
NGC1808 Saskia et al. 1990), there are few galaxies which 
have been systematically surveyed for complete samples 
of SNR. Of these, the best candidates are the Magellanic 
Clouds (eg. Matthewson et al. 1983) and M33 (eg. Long et 
al. 1990). This paper reports on our initial study of a sample 
from this latter set. 



2 THE SAMPLE AND THE OBSERVATIONS 

Our sample of objects was chosen from the confirmed SNR in 
the catalogue of Long et al. (1990). Since Oliva et al. found 
a correlation between the 1.644/tm [Fell] line and H/?, we 
imposed a cut using the Long et al. Ha surface brightnesses 
of 5xl0" 19 Wm" 2 arcsec" 2 . The limit in 1.257/tm [Fell] sur- 
face brightness is then ~10% of this value (using the rela- 
tionship between 1.644/tm [Fell] and H/?, the fact that the 
1.257/tm/1.644 / (tm ratio is ~ 1.35, and a standard case B 
recombination Ha/H/? ratio, though we note this latter may 
not be strictly valid). A second requirement was that the 
sources should be relatively compact given that we were ob- 
taining spectroscopic data. We imposed a limiting diameter 
of ~10arcsec (~35pc in table 3 of Long et al). The sample 
discussed in this paper is incomplete: the complete surface 
brightness limited sample derived from Long et al. is given 
in Table 1, and the ten SNR we actually observed are noted, 
together with the date of observation and integration time. 
Ideally for our purposes, we would have used a radio se- 
lected sample, chosen such that it spanned a wide range in 
radio luminosity. However, at the time of our observations 
such a catalogue was unavailable, though one has now been 
published (Duric et al. 1993). In practice, the sources we 
observed do span a reasonable range of radio luminosity as 
will be shown later. 

All objects were observed with the common-user array 
spectrometer CGS4 on UKIRT (Mountain et al. 1990) on 
the nights of 13, 14 and 15 November 1992. Conditions were 
assessed as being photometric from the measurements of the 
flux standards described below. CGS4 utilises a SBRC InSb 
62x58 array, sensitive from 1-5/tm. The spatial resolution 
was 3" per pixel, and the slit-width was also 3". In order to 
obtain fully sampled data, the array is stepped 4 times over 2 
pixels: this also enables single bad pixels to be removed from 
the data (Puxley, Beard & Ramsay 1992). The central wave- 
length for all observations was 1.26/tm. We chose to observe 
the 1.257/tm a 6 D 9 / 2 -a 4 D 7 / 2 [Fell] line since it is inherently 
stronger by a factor of 1.35 (Nussbaumer & Storey 1988) 
than the 1.644/tm a 4 F 9 / 2 -a 4 D7/ 2 line, allows simultaneous 
observation of Pa/?, is significantly less affected by night sky 
emission bands than the H window (though a strong and 
variable oxygen airglow line is centred at 1.27/tm - eg Ram- 
say, Mountain & Geballe 1992) and extinction is not a prob- 
lem given the inclination of 55° of M33. The objects were 



observed in 2nd order, with a resolution of approximately 
400. 

As atmospheric standards we used the late-type A stars 
SA054655 (AO), BS553 (A5V) and HD3029 (A3), from 
which the Pa/? absorption feature was removed by fitting 
a Lorenzian absorption profile. All standards were observed 
at similar air mass to the target SNR, eliminating the need 
to make airmass correction in our flux calibration. We ob- 
served one of these standards once every 1-2 hours. BS553 
and HD3029 are also flux standards in the UKIRT standards 
list, with magnitudes of J = 2.50 and J = 7.263 respectively. 
Using these values we obtain J = 5.17 ±0.05 for SA054655. 
Wavelength calibration was by means of krypton and xenon 
arc lamps. Comparison of the fluxes of the standard stars 
from repeated observations over the three nights provides 
a measure of the accuracy of the overall photometry. From 
this we estimate that the photometric errors on any given 
single spectrum are <10%. 

The slit was oriented north-south for all observations. 
We pointed at the positions given by Long et al., though 
for DDB-2 we also searched near the location given by 
D'Odorico, Dopita & Benvenuti (1980). In all cases in which 
we detected emission, the peak was found to be within ap- 
proximately 1-2 arcsec of the quoted positions. The objects 
observed are discussed more fully below. For all but DDB- 
2 and 013110+30210 we are confident that we observed at 
the location of the optical SNR. For those objects that were 
known from the optical data to be extended, we stepped the 
telescope a full slit width across the object until no flux was 
detectable in the same integration time as that for the cen- 
tral position (generally 320 or 640 seconds - see Table 1). 
We did not map further in those cases where we failed to 
detect [Fell] emission in the central position. In most cases 
we also integrated on the central position for a further con- 
siderably longer time in order to obtain good signal-to-noise 
on the less luminous SNR and to detect further weaker lines 
in the brighter SNR. These repeat observations also allowed 
us to check the accuracy of our photometry. For the com- 
pact sources the differences are less than 30%, for the more 
extended ones it is closer to 100%. We conclude therefore 
that our overall photometry for all sources is accurate to 
this level. Future surveys of these objects seeking accurate 
photometry should ideally use narrow band imaging tech- 
niques. 

The reduction procedure for these data followed a stan- 
dard procedure. First, the sky frames were subtracted from 
the corresponding object frame. It was found in practice 
that up to 10% of the original sky line flux was left after 
this subtraction. Therefore, a polynomial fit to the residual 
sky emission lines in this object-sky pair was then made, 
and the resulting fit removed to correct for this. The indi- 
vidual rows containing the positive and negative beams were 
then extracted. The data were then coadded, and corrected 
for atmospheric transmission by division by a standard star. 
Lastly spectra were then corrected for the assumed black- 
body shape of the standard. Effective temperatures for the 
standards were assumed to be 8200K, 8720K and 9520K for 
BS553, HD3029 and SA054655 respectively. 

No extinction correction has been applied to this data. 
From Blair & Kirshner (1985), it can be seen that the extinc- 
tion at J is likely to be small and hence is entirely negligible 
for the purposes of calculating line ratios. For absolute line 
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luminosities, the correction is mostly <10%. The observed 
fluxes are given in Table 2 and the spectra of those SNR 
with positive detections are shown in Figure f. All spectra 
shown are of the single long exposures on the peak positions 
only. The error bars shown represent both overall photo- 
metric error on the spectrum, but also provide an indica- 
tion of those regions where the sky transmission is less, or 
where the residual contamination from sky emission lines is 
greater. Since the positions observed away from the peak of 
the remnants were often of low signal-to-noise, weaker lines 
were not generally detected in these locations. Generally the 
line strengths listed in Table 2 will reflect integrated inten- 
sities for the strong lines but only peak intensities for the 
weaker ones, and comparison of ratios of these is not strictly 
valid. Therefore, the comparison of line ratios in what fol- 
lows (especially section 3.2) applies to lines measured in the 
peak position alone. 

3 RESULTS 

3.1 Comments on Individual Remnants 

DDB-2: The centroid of this object as given by Long et al. is 
~5arcsec east of the coordinates given by Duric et al. for the 
radio centroid. D'Odorico et al. (f 980) give a position some 
23arcsec south of the Long et al. one (this is the largest 
discrepancy between these two references for objects that 
are common to both that we actually observed). We tried 
both the Long et al. and D'Odorico et al. positions with- 
out detecting any [Fell] emission. However, since we did not 
map beyond these two central positions it is possible that 
we narrowly missed the object when observing the Long et 
al. centroid (the object itself is compact with a diameter of 
<4arcsec). 

DDB-4: Weak [Fell] emission was detected on the central po- 
sition. No emission was evident either 3arcsec east or west. 
There is a very small offset between radio and optical po- 
sition well within the claimed accuracy of the Long et al. 
coordinates (2arcsec). Pa/? was detected at ~ 50% of the 
1.257/tm [Fell] flux. 

013022+30244: No emission was evident. Again, the cen- 
troid given by Long et al. is some 5arcsec east of the radio 
centroid. The 20cm radio map of Duric et al. also shows con- 
siderable emission south-west of this location, though only 
very diffuse nebular emission appears to be present there. 
Again, it is possible that we may have missed this very com- 
pact optical source. 

DDB-6: Very small offset between radio and optical position 
of 2arcsec, though in the same sense as above (radio posi- 
tion west of optical). [Fell] emission is extended: we found 
detectable emission in 3 adjacent slit positions. Very weak 
Pa/? emission was also evident. 

DDB-7: Radio and optical positions agree within the errors. 
Faint emission evident beyond the central slit position, but 
bright core of object is small. Weak Pa/? emission evident. 
DDB-8: Radio and optical positions agree within the errors, 
though we find peak of [Fell] emission to be in the slit po- 
sition 3arcsec west of this. Emission is more diffuse than in 
the case of DDB-7. Pa/? in this source is stronger at ~30% 
of [Fell] on the core, and brighter than the [Fell] line in the 
west. This HI emission may be background. 



DDB-9: Again we find [Fell] peaks to the west of the ra- 
dio/optical centroids. Long et al. note that the object ap- 
pears circular with the western edge being brighter. Moder- 
ately strong Pa/? emission across the SNR. 
013109+30225: Not detected in the radio survey of Duric 
et al. Faintest detection in our survey. No evidence for ex- 
tended emission along the slit on the central position, so 
only one map position was observed. 

013110+30210: Again no detected radio emission. No de- 
tectable [Fell] emission though we do detect the star seen to 
the south of the optical source. We are therefore confident 
that our positioning for this object is correct. 
DDB-11: We find the peak to lie to the west of the ra- 
dio/optical centroid. Long et al. note that the [OIII] emission 
also peaks north-west of the centroid. Weak Pa/? to the east 
of the core, none on the brightest [Fell] spot. 

For the three non-detections we can be absolutely con- 
fident of our positioning for only one (013110+30210). It is 
also worth noting that the compactness of the [Fell] emission 
does not necessarily follow the radio, or the quoted optical 
diameters from Long et al. We shall return to this issue in 
section 3.4. 



3.2 Physical Conditions Within the SNR 

We can use the ratios of the observed [Fell] lines to place 
constraints on the physical conditions present within the Fe 4 " 
emitting zone. The atomic data of Nussbaumer & Storey 
(1988) and Pradhan & Zhang (1993) were used to interpret 
these ratios. The much larger collision strengths reported by 
Pradhan & Zhang compared to Nussbaumer & Storey (1980) 
reflects the fact that their computations take account of the 
many resonances that exist in the excitation of Fe 4 " more ac- 
curately. Since the method employed should in principle be 
more accurate than that used in obtaining the older Nuss- 
baumer & Storey (1980) data we will use their data in this 
analysis. 

The 1.257/tm line luminosities of the brighter remnants 
are > I0 29 W . For typical densities and temperatures, the 
luminosity of all the low lying [Fell] lines together is about 
ten times larger than this. Since a supernovae injects around 
I0 44 J into the interstellar medium, if the typical [Fell] emit- 
ting time is only of the order of 1 4 years (see below) then 
it is easy to see that Fe 4 " alone is not an important coolant. 

All of the [Fell] lines detected by us in these observa- 
tions arise within the a 4 D-a 6 D multiplet (see Nussbaumer 
& Storey 1988 for a complete list of the lines in this mul- 
tiplet). All of our detections are consistent within the er- 
rors with the available atomic data, as are the upper lim- 
its on other lines in this mulitplet, with the exception of 
the 1.279/tm a 4 D 3 /2-a 6 D 3 / 2 transition. This line is blended 
with Pa/? in many cases, and also suffers from being coinci- 
dent with strong night sky lines of OH and O2. We therefore 
attribute the discrepancy in derived electron density using 
this line to our inability to both fit the line accurately and to 
residual sky contamination. Therefore, all electron density 
determinations are derived using the ratio of the 1.294/tm 
a 4 D 5 /2-a 6 D 5 /2 to 1.257/tm a 4 D7/ 2 -a 6 D 9 / 2 transitions. Lim- 
its are derived by measuring the flux in a 2 pixel wide seg- 
ment at the correct wavelength for the 1.294/tm line. 
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Both Oliva, Moorwood & Danziger 1990 (for the SNR 
RCW 103) and Rudy et al. 1992 (for the emission line star 
MWC 922), reach substantially the same conclusions on the 
comparison of observations and theory using the older col- 
lision strengths of Nussbaumer & Storey (1980). Since the 
new strengths of Pradhan & Zhang are substantially larger 
in most cases, the similarity of these conclusions indicates 
that most of the changes are due to an overall scaling of 
all the collision strengths (see the discussion in Oliva et al. 
1990). To see the effect the changes have had, we consider 
also the electron densities as derived from the optical data of 
Blair & Kirshner (1985) and Smith et al. (1993). The [SII] 
transitions at 6717A and 673lA are sensitive to densities 
in the range n e = 100 — 10 4 cm -3 . Table 3 shows the den- 
sities derived using this method compared to the densities 
derived using the [Fell] ratios. The error in the derived n e 
comes from the effect of varying the electron temperature (a 
negligible effect for the [Fell] densities). Since we have ob- 
served the brightest objects in the Long et al. catalogue, the 
observed signal-to-noise in the optical lines is always large. 

From these results we can see that the derived elec- 
tron densities in general are less from the [SII] ratios than 
from the [Fell] ratios. This fact was first noted by Oliva, 
Moorwood & Danziger (1989) for a small sample of galactic 
and LMC SNR using the older Nussbaumer & Storey (1980) 
atomic data. However, we see a smaller offset than they did. 
We therefore also checked the results of Oliva et al. using the 
newer atomic data for Fe 4 " as well. They used the density 
sensitive ratio 1.599/tm a 4 D 3 / 2 -a 4 F7/ 2 to 1.644/tm a 4 D7/ 2 - 
a F 9 12 ratio to derive electron densities in a small sample of 
SNR in our galaxy and the Large Magellanic Cloud. Using 
the Pradhan & Zhang collision cross-sections, the predicted 
values of n e are all less than given by Oliva et al. These val- 
ues are also given in Table 3. This reduces the discrepancy 
between the inferred electron densities from the [Fell] and 
[SII] ratios they see, though it does not eliminate it, just as 
is the case for our observations of the M33 SNR. Given the 
change that has occurred in moving from the Nussbaumer & 
Storey (1980) data to the Pradhan & Zhang data we must 
caution that the remaining discrepancies may also still be 
due to inaccuracies in the collision cross sections. However, 
their rationale for this discrepancy, the low critical density 
for the [SII] lines (and in particular the 6717A line), is also 
perfectly valid. It would be extremely valuable to obtain 
further data on IR transitions of [Fell] in bright galactic 
remnants with the current generation of IR spectrometers 
to further test the reliability of the atomic data. 

Lastly, we note that there is a relatively wide span in 
SNR [Fell] luminosity in our sample, but little evidence for 
there being a wide range in density. We therefore conclude 
that the ambient density plays only a minor role in deter- 
mining the [Fell] luminosity of the SNR in our sample. 



3.3 Comparison of [Fell] and other observed 
emission 

Both Smith et al. (1993) and Blair & Kirshner (1985) have 
obtained optical spectra of the SNR in M33 in the wave- 
length range ~ 4000A-7500A. It is worthwhile comparing 
the fluxes arising from different species, to seek possible cor- 
relations. It has already been noted (for example Smith et 



al.) that the forbidden line fluxes tend to be correlated if 
they are not strongly dependent on the physical conditions 
in the SNR. Unfortunately, Blair & Kirshner observed only 
some of those objects in the D'Odorico et al. sample. We 
therefore compare our data with that of Smith et al. who 
observed all the SNR under consideration here. Smith et al. 
only observed at one slit position per object, but our own 
data suggests that the sources are sufficiently compact that 
the bulk of the emission does lie in a narrow slit. Smith et 
al. also aligned their slit over the brightest regions of the 
SNR: the line fluxes as measured by them for the objects in 
our sample are given in Table 4. Also given there are opti- 
cal sizes from Long et al. (1990), radio sizes and fluxes from 
Duric et al. (1993) and velocity information from from Blair, 
Chu & Kennicutt 1988 and Blair & Davidsen 1993. 

In Figures 2a-d we show how the 1.257/tm [Fell] flux 
correlates with 6300A [OI], 673lA [SII], 6365A Ha fluxes 
and 6cm radio flux, and present linear correlation coeffi- 
cients for all these tests in Table 5. In this analysis, objects 
with only upper limits have been excluded. The 673lA [SII] 
line was chosen since it has a much higher critical density 
than the 6717A line. Clearly there is a well defined correla- 
tion between [Fell] emission and the optical lines, but only 
a weak correlation between the forbidden line fluxes and 
the measured radio continuum. The discrepancies that exist 
between our [Fell] data and the optical emission line data 
are probably due largely to the different effective apertures 
being used. 

Most of these correlations (or lack of them) can be un- 
derstood in terms of the source of the emission. Generally, 
the optical and infrared line emission in SNR arises in ra- 
diative shocks, whereas the radio emission arises from syn- 
chrotron radiation. In particular, the fluxes of [SII], [Fell] 
and [OI] should all be strongly correlated since they all arise 
in a region in which hydrogen is largely neutral. This is true 
both of shock excited emission, and of photoionisation (as 
is seen in the Crab Nebula for example - e.g. Fesen & Kir- 
shner 1982, Henry & MacAlpine 1982). By contrast, [Oil], 
with a similar ionisation potential to hydrogen, comes from 
a zone in which hydrogen is at least partially ionised. In 
terms of shock models, the various regions that we can ob- 
serve emission from are the immediate post-shock region 
(where hydrogen is fully ionised, [OIII] emission arises and 
iron if present is in a higher ionisation state), a transition 
region in which gas recombines rapidly (the [Oil] region) 
and a larger volume of gas in which ionised stages of iron 
are rapidly converted into Fe 4 " by charge exchange reactions 
with the now dominant neutral hydrogen (see, for example, 
Oliva et al., Van der Werf et al. 1993). 

The one extremely discrepant point in Figure 2d, DDB- 
9, may have a different nature from the other remnants, since 
it is also a bright x-ray source (eg Long et al. 1990: DDB-7 
has also been detected in x-rays but shows no signs of being 
different from the other sources in other regards) However, 
from our current data, we cannot say anything further. 



3.4 Correlation of [Fell] luminosity and SNR 
properties 

In trying to determine how the [Fell] luminosity correlates 
with the age of the remnant, we must first decide on a 
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suitable observable quantity that defines this latter param- 
eter. In principle, virtually all of the remnants in our sam- 
ple should be either in an adiabatic or radiative expansion 
phase. Long et al. (1990) found that N "(< D) oc D 2 1 (num- 
ber of SNR with diameter less than Dpc) for M33, which 
should be compared with an exponent of 2.5 for adiabatic 
expansion and I for free expansion. Hence, we will assume 
that the standard Sedov solution can be applied to the evo- 
lution of these remnants (see, eg. Lozinskaya 1992 for a fuller 
description of the evolution of SNR). Although this is un- 
likely to be accurate at early times, it should provide a rea- 
sonable approximation at these late epochs. We have used 
this solution to derive approximate ages, in the first instance 
by using the relationship between age and remnant diame- 
ter, but also by considering the relationship between shock 
velocity and age as well. According to the Sedov solution, 
these quantities are then related by 

R s oc (E /n ) 2 t 0A v s = 0AR s /t, 

where R s and v s are the radius and velocity of the shock, t 
is the remnant age, Eo the initial energy of the supernova 
and no the density of the amibient medium. For typical pa- 
rameters (Eo = I0 44 J, n = Icm" 3 ), R s ~ 0.34i° r 4 pc. 

Although it is likely that the remnants as a whole 
are still undergoing adiabatic expansion, radiative shocks 
are prominent in these objects. One possibility is that the 
[Fell] lines reach maximum luminosity when the remnant as 
a whole undergoes the transition from adiabatic to radia- 
tive expansion. Detailed models (eg. Cioffi & McKee 1988) 
show this to be a relatively slow process, and indeed pre- 
dict the SNR luminosity to peak just as the radiative phase 
is reached. Typically, this would occur after ~ I0 5 years in 
their models, rather older than the brightest remnants that 
we see. 

There are many different estimates of the sizes of these 
SNR. Long et al. give optical diameters based on the ge- 
ometric mean of the dimensions of the SNR found in their 
optical images. Duric et al. quote diameters based on the de- 
convolution of their radio data. These do not always agree 
as shown by Table 4. Blair & Davidsen (1993) find different 
results to either of these groups based on analysis of Hubble 
Space telescope imaging data of DDB-2 and DDB-8. Blair, 
Chu & Kennicutt (1988) also give different diameters as de- 
termined from their echelle spectroscopy. Similarly, higher 
resolution VLA data (Goss & Viallefond 1985) shows that 
DDB-7 is substantially smaller than either Long et al. or 
Duric et al. report. We also find a smaller diameter for this 
source (~ 3.5") from an [Fell] image obtained with IRCAM3 
at UKIRT. 

Given these problems with determining the 'true' diam- 
eters of these SNR, it is not altogether surprising that the 
relation between [Fell] flux and optical diameter shown in 
Figure 3a is so poor (using radio diameters is no better as 
the correlation coefficients in Table 5 show). We therefore 
rule out using diameter as an accurate measure of age. The 
better correlation found between the optical emission lines 
and the optical radii is effectively a selection effect. Smith et 
al. (1993) align their slit to 'include the brightest portions 
of the objects' from the original Long et al. images. It is not 
surprising then that the measured flux correlates with the 
original data, but it does not discount the possibility raised 



by Blair & Davidsen (1993) that the measured diameters of 
all these sources are wrong. 

By comparison, in Figure 3b we show the correlation 
between [Fell] flux and velocity for the remnants. The cor- 
relation here is much stronger (Table 5), with the expected 
trend that the most luminous sources have the broadest Ha 
lines. For typical values of the parameters (initial energy, 
density of the interstellar medium), the Sedov solution (eg. 
Lozinskaya 1992) implies that the brightest remnants in this 
sample are ~ I0 4 years old. We also note the poor correla- 
tion between optical radius and velocity width (Figure 3c): 
the measured values of the diameters and velocity widths 
for the two most [Fell] luminous remnants in our survey are 
completely incompatible with the Sedov solution. These dis- 
crepancies are not the result of obvious density differences, 
since the data from the [SII] lines shows that most of these 
SNR have similar densities. Hence our conclusion on the re- 
liability of the measured diameters as an age probe is borne 
out by this comparison as well. Unfortunately, since we do 
not have measured velocity widths for any of the most com- 
pact optical remnants, we cannot tell if they are also gen- 
uinely young using this measure of age. 

Therefore, from this dataset alone, we cannot rule out 
either of the two possibilities: first, that the forbidden line lu- 
minosity is largest for the youngest remnants with the most 
energetic shocks; secondly, that the luminosity peaks at the 
point which most likely represents the slow transition from 
adiabatic to radiative expansion. We can say from the ob- 
served weakness of [Fell] emission in many of the younger 
remnants in our galaxy and the LMC, that there must be 
some lower cut-off as well as a higher one (Hester et al. 1990 
find a 1.644/tm luminosity for the Crab remnant ~ I0 27 W, 
and we similarly find little evidence of [Fell] emission in the 
youngest of the LMC remnants we have studied). Hence the 
idea that [Fell] luminosity may peak sharply with the age 
of the remnant is consistent with our results but further 
observational data is required to confirm it. 



4 DISCUSSION 

From our incomplete sample of SNR in M33, we have demon- 
strated the fundamental similarities between the infrared 
[Fell] line strengths and those of well studied optical tran- 
sitions with equivalent ionisation potential such as [SII]. 
Therefore, [Fell] lines are good tracers of SNR properties. 
This is especially true for regions of high extinction (given 
their longer wavelength) and where HII regions coexist in 
the beam with SNR (since [Fell] emission from HII regions 
is essentially negligible). 

It is useful to compare the typical luminosities that we 
find here with the values found for the 1.644/tm line by other 
groups. The only well studied older remnants are in our own 
galaxy and the LMC (see Oliva et al. and references therein). 
Typical bright remnants there have equivalent 1.257/tm lu- 
minosities of ~ I0 30 W, in reasonable agreement with the 
values we find for M33. These peak luminosities also appear 
to arise from remnants which are typically a few thousand 
years old. In a future paper we will present our own findings 
on the [Fell] luminosities and structures of a larger sample of 
LMC remnants. It seems likely from the work completed to 
date that the dominant factor in determining the luminosi- 
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ties for these evolved remnants is their age. It is unlikely that 
metallicity is the cause (both the LMC and M33 have simi- 
lar metallcities, yet the brightest [Fell] source in the LMC, 
N49, is almost ten times brighter than DDB-8). Similarly, 
variations are seen in the density of the bright remnants 
(as is the case here for DDB-7 and DDB-8) so again this is 
unlikely to be a major factor. The strong correlation that 
we see between luminosity and bulk velocity indicates that 
age is the dominant mechanism within a galaxy, and this is 
likely to be true when comparing different systems as well. 

Given these conclusions, it is worth discussing briefly 
the nature of the compact radio sources seen in many 
nearby starburst galaxies, and commonly assumed to be 
SNR (though see Ulvestaad & Antonucci 1994). The com- 
pact sizes of these 'remnants' naturally leads to the assump- 
tion of a young age for them. However, there are two prob- 
lems with this simple interpretation. First, the narrow width 
of the [Fell] emission detected by Greenhouse et al. in M82 
is rather puzzling if this is true. Secondly, the radio data 
shows no strong evidence for a rapid fading of these sources 
as is seen in isolated radio supernovae (eg. Ulvestaad & An- 
tonucci 1994, Huang et al. 1994). If we assume that these 
SNR are young, then we must further conclude that the 
emission mechanism in these remnants is a different process 
to the radiative cooling of filaments that we see in these 
older remnants. One possibility has received considerable 
attention to model the luminous radio emission seen from 
some extragalactic SNR. In this model the supernova blast 
wave overtakes a dense circumstellar shell, and it is largely 
the excitation of this shell that we see at early times (eg. 
Chevalier & Frannson 1994). In the high pressure environ- 
ment in the nuclei of starburst galaxies, it seems likely that 
such shells would be confined even closer to the progenitor 
star, and hence may lead to even greater emission. Even 
this model has some difficulties however, since it predicts 
broad emission lines should be present. Further study of the 
properties of remnants in a high pressure environment, both 
theoretical and observational, is highly desirable. 
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Table 1. The complete diameter and surface brightness limited sample of SNR in M33. The Ha surface brightness is given in units of 
10- 15 Wm" 2 arcsec 2 . Diameters and surface brightness values are taken from Long et al. (1990) 
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Table 2. Observed line fluxes seen in the M33 SNR. No extinction corrections have been applied. Lines marked with an * were only 
detected at the peak of the emission, as discussed in the text. All fluxes are in units of 10- 18 Wm _2 . 



Object n e 
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Table 3. Derived electron densities from [SII] (Smith et al. 1993) and [Fell] transitions. The error bars for the [SII] data reflect the 
uncertainty in the electron temperature rather than any observational error. The observed signal to noise is large enough for the [SII] lines 
that this is valid. The variation with electron temperature for the [Fell] lines is negligible for most values of n e . The electron densities 
for the Oliva et al. sample derived from the [SII] lines are taken from that source. The densities derived from the [Fell] lines have been 
calculated as described in the text. 
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Table 4. Published values for line fluxes, sizes, velocity data and radio fluxes taken from Smith et al. 1993, Long et al. 1990, Duric et 
al. 1993, Blair, Chu &l Kennicutt 1988 and Blair &l Davidsen 1993 as noted in the text. Limits are given where appropriate, and a blank 
is left where no published data exists. 
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Table 5. Linear correlation coefficients for the various measured parameters of the SNR discussed in the text. 
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Figure 1. J-band spectra for DDB-4, DDB-6, DDB-7, DDB-8, DDB-9, DDB-11, 013109+30225. Spectra are presented for the peak 
position only. 
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Figure 1. 
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Figure 2. Scatter plots for 1.257/xm [Fell] against (a) 6300A [OI] (b) 673lA [SII] (c) Ha and (d) 6cm radio flux density. Upper limits 
are shown as arrows, and are 3<r limits. 



12 S.L. Lumsden and P.J. Puxley 




5 10 

Optical diameter (arcsec) 




200 400 600 800 
Bulk Velocity Width (krns -1 ) 



o 

o h 
co 



n — r 



o l± i i i_ 





_i i i i i =i 

5 10 



Optical diameter (arcsec) 



Figure 3. Scatter plots for 1.257/im [Fell] against (a) optical diameter and (b) bulk velocity width (see Blair et al. 1988 for a discussion 
of the velocity widths). Also shown is the relationship between velocity width and radius. 



